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EFFICIENT CHANNEL ESTIMATION IN A 
DIGITAL COMMUNICATIONS SYSTEM 



BACKGROUND OF THE INVENTION 



1. Field of the Invention 



[0001] 



The present invention relates generally to the field of digital communications; 



and, more particularly, to an efficient channel estimation method and apparatus for a 
5 GSM/EDGE digital communications system which utiUzes special properties of GSM/EDGE 
training sequences and leads also to a practical joint optimization of synchronization and 
equaUzer window sizing. 



10 systems, including those systems which operate in accordance with the Global System for 
Mobile Communications (GSM) specifications. ISI is caused by time dispersion in the 
transmission channel over which a signal is transmitted, and adversely affects the quality of the 
received signal. In effect, ISI causes distortion of the transmitted signal which, in turn, causes 
symbol errors in the received signal such that it becomes difficult for the receiver to determine 

1 5 what data was actually sent. 
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2. Description of the Prior Art 



[0002] 



Inter-symbol interference (ISI) is an important problem in digital communications 
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[0003] As is well-known in the art, the usual way to compensate for ISI in a GSM system 

is to provide a channel estimation-based equalizer in the receiver. Basically, a model or estimate 
of the propagation channel over which a received signal was transmitted is created; and the 
equalizer then uses that information to estimate the sending symbols that best fit the received 
5 signal 

[0004] EDGE (Enhanced Data rates for Global Evolution) is an interface mode which has 

recently been developed for GSM Networks. EDGE'S principal features include new modulation 
and coding schemes which increase data capacity and speed in the air interface. EDGE is fully 
based on GSM and uses the same TDMA (Time Division Multiple Access) frame structure as 
10 GSM, such that it allows GSM operators to use existing GSM radio bands to offer wireless 
multimedia-based services and applications. 

[0005] In GSM/EDGE systems, the performance of the equalizer in combating ISI 

depends heavily on the quality of the channel estimation; and the quaUty of the channel 
estimation is, in turn, highly sensitive to the acciu-acy of burst synchronization (the term 
15 "synchronization" in this document signifies "burst synchronization") and the size of the 
equalizer window. 

[0006] To effectively combat ISI, the span of the equaUzer window must be large enough 

to cover the maximum delay spread of the channel. However, an over-sized equalizer window 
will result in performance losses in channels with shorter delay spread due to inaccurate channel 
20 estimation, inadequate synchronization and increased noise contribution. In addition, 
synchronization can be optimized only with the knowledge of the equalizer window to capture 
maximum energy in the window span. In general, the interdependency of synchronization and 
equalizer window sizing makes efficient optimization difficult to achieve under different channel 
conditions. Joint optimization was deemed as too expensive in practical implementation. 

2 
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[0007] Because of this interdependency, a careful compromise has to be made for a 

traditional equalizer with fixed window span so that the size of the window can provide adequate 
performance in long delay spread channels, such as Hilly Terrain (HT), without losing too much 
in short delay spread channels, such as Typical Urban (TU). 
5 [0008] When the equalizer window size is set, there are two approaches for burst 

synchronization, a correlation-based approach and a least square error (LSE) approach. As will 
be described more fully hereinafter, in current GSM receivers, synchronization is done by a 
correlation-based approach in which the central 16 symbols of a known 26-symbol training 
sequence is correlated with the central 16 symbols of the training sequence in the received signal, 

10 The synchronization is determined by either the centers of gravity of the correlation or the 
maximum correlation energy in the equalizer window. In the LSE approach, for each possible 
synchronization point, an estimation of the chaimel is made, and synchronization is determined 
by comparing the expected and received training sequence with least square error criterion. 
[0009] Neither of these approaches, however, is fully satisfactory. The correlation-based 

15 algorithm suffers from performance degradation due to inaccurate synchronization, especially in 
long dispersive channels, while a straightforward LSE-based algorithm suffers from a high 
degree of computational complexity (mainly due to multiple channel estimation). 



SUMMARY OF THE INVENTION 
20 [0010] The present invention provides an initial channel estimation method and apparatus 

in a GSM/EDGE digital communication systems which utilizes special properties of 
GSM/EDGE training sequences to achieve a more efiicient channel estimation. A two- 
dimensional recwsion is further estabhshed to enable practical joint LSE optimization of burst 
synchronization and adaptive equalizer window sizing, 

3 
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[0011] More particularly, an embodiment of the present invention provides a method for 

efficient estimation of the transmission channel in a digital communications system which 
operates in accordance with GSM/EDGE specifications. The system has a receiver which 
receives a signal transmitted over the transmission channel and which includes a channel 
5 estimation-based equalizer that uses an estimate of the transmission channel to equahze the 
received signal. A received signal burst includes a 26-symbol training sequence. The method 
includes using any consecutive 16-symbol segment of the 26-symbol training sequence to 
estimate, without matrix operation, the transmission channel at different possible synchronization 
points for different equalizer window sizes. 

10 [0012] The present invention provides a more efficient channel estimation method by 

recognizing and utilizing two previously unexploited properties of GSM/EDGE training 
sequences which arise from the cyclic prefix structure of the training sequences. In particular, the 
present invention recognizes that the 26-symbol GSM/EDGE training sequences are, within 
certain ranges, both shift invariant and order invariant. The property of being shift invariant 

15 enables a channel estimation to be carried out with delayed (shifted) training sequence segments; 
which, in turn, permits: L The ISI corrupted leading training sequence symbols to be avoided as 
much as possible in long dispersive channels; and 11. The leading taps of a channel to be 
estimated using the same training sequence segments, regardless of the size of the equalizer 
window. The property of being order invariant permits channels with different time dispersion, 

20 from 1-8 symbols (as complex polynomials of 1-8*^ order) to be estimated without matrix 
inversion if any consecutive 16-symbol segment of the training sequence is used. This results in 
a significant reduction in computational complexity in channel estimation for varying window 
size which can be adaptive to the actual delay spread of the channel. Combining these two 
properties, a two-dimensional recursive relation can be established for different synchronization 

4 
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points and tap positions, which makes a practical LSE joint optimization of synchronization and 
equalizer window sizing possible. 

[0013] In general, the present invention provides an efficient channel estimation for 

GSM/EDGE systems which, beside direct utilization, is fundamentally important to achieving a 
5 practical implementation of joint LSE synchronization and equalizer window sizing. 

[0014] Yet further advantages and specific details of the present invention will become 

apparent hereinafter in conjunction with the following detailed description of presently preferred 
embodiments. 

10 BRIEF DESCRIPTION OF THE DRAWINGS 

[0015] Fig. 1 schematically illustrates the structure of a TDMA frame used in a digital 

communications system according to the GSM/EDGE standard; 

[0016] Fig. 2 schematically illustrates a cyclic prefix structure of GSM/EDGE training 

sequences; 

1 5 [0017] Fig. 3 is a block diagram of a simplified baseband model for a transmitter-receiver 

chain to illustrate data recovery signal processing in GSM/EDGE digital communications 
systems; 

[0018] Fig. 4 is a flow chart illustrating channel equalizer processing in a GSM/EDGE 

digital communications system; and 
20 [0019] Fig. 5 illustrates a two-dimensional recursive relation in channel estimation for 

different synchronization points and channel tap positions according to a presently preferred 
embodiment of the invention. 
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DETAILED DESCRIPTION OF PRESENTLY PREFERRED EMBODIMENTS 

[0020] Fig. 1 schematically illustrates the structure of a TDMA frame 20 used in a 

mobile telecommunications system according to the GSM standard. Frame 20 has a duration of 

4.615 ms and accommodates 8 information channels (time slots) 21, designated 0-7 in Fig. 1. 

5 Each of the 8 time slots has a duration of 0.577 ms and contains a 148-bit signal portion 30 and 
an 8.25-bit guard portion (not shown) which functions to maintain separation between signals in 
adjacent time slots. The 148-bit signal portion is generally referred to as a normal burst and 
comprises a first 3-bit tail bit section 22, a first 57-bit coded data section 23, a first 1-bit 
housekeeping bit section 24, a 26-bit training sequence section 25, a second 1-bit housekeeping 

10 bit section 26, a second 57-bit coded data section 27 and a second 3-bit tail bit section 28. 

[0021] The 26-bits comprising the training sequence section 25 in a GSM digital 

communications system, is typically viewed as being divided into a centrally located portion 31 
of 16-bits (sometimes referred to as the midamble portion), and side portions 32 of 5-bits each. 
Alternatively, in a GSM/EDGE digital communications system, the 26-bit training sequence 25 

15 can be also viewed as including a 16-bit white sequence 37 and a 10-bit cyclic prefix 38 as is 
illustrated in Fig. 2. 

[0022] As mentioned previously, inter-symbol interference can occur in digital 

communications systems as a result of time dispersion in the transmission channel over which a 
signal is transmitted; and in GSM/EDGE systems, the usual way to compensate for ISI is to 
20 provide a channel estimation based equalizer in the receiver. 

[0023] Fig. 3 is a block diagram of a simpUfied baseband model of a transmitter-receiver 

chain, generally designated by reference number 40, to illustrate data recovery signal processing 
in a GSM/EDGE digital communications system. In Fig. 3, the transmitter side is shown at 41 
and the receiver side is illustrated at 42. A signal is transmitted from the transmitter side 41 to 
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the receiver side 42 over a transmission channel 43. The transmission characteristics of the 
transmission channel frequently change due to a variety of factors, including the motion of the 
mobile terminals, the fluctuation of multi-path (time dispersive) propagation channels and 
variant interferences, which are generally designated in Fig. 3 by the noise component 44 
5 introduced into the transmission channel As is well-known to those skilled in the art, the multi- 
path channel and the noise component can adversely affect the quality of the received signal, 
e.g., cause intersymbol interference; and this necessitates that the received signal be corrected, 
typically by means of a channel estimation-based equalizer. 

[0024] Specifically, as shown in Fig. 3, at the receiver side, the received signal, which 

10 includes the transmitted signal with intersymbol interference and the noise component, is, after 
being filtered by filter 45, sent to an equalizer 46 which attempts to correct the received signal 
utilizing an estimate of the transmission channel as shown at 47. The chaimel estimation, in tum, 
is a function of burst synchronization and equalizer window sizing as shown at 48. The 
performance of the equalizer can be improved significantly if, in an advanced receiver, the 
15 equalizer window size is adaptive to the actual delay spread of the channel Data representing an 
estimate of the actual transmitted symbols, as determined by the equalizer, is then output for 
fiirther processing by the receiver. 

[0025] Fig. 4 is a flow chart which illustrates, in somewhat greater detail, the steps of the 

channel equalization processing which takes place in an advanced receiver 42 of Fig. 3. 
20 Specifically, after filtering the received signal at step 50, synchronization and equalizer window 
sizing according to the delay spread of the actual channel is performed as shown at step 52. An 
estimate of the transmission channel is then made at step 54, and a hypothetical signal provided 
in step 60 is applied to the estimated channel at step 56. An error comparison is made between 
the hypothetical signal applied to the estimated channel from step 56 and the synchronized 

7 
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received signal at step 58. If the result of the error comparison is not at a minimum value (not 
min output of step 58) a revised signal hypothesis is made in step 60, the new revised signal is 
applied to the hypothetical channel at step 56, and a new comparison is made at step 58. The 
comparison process is continued until the result of the comparison is at a minimum value (min 
5 output of step 58), at which time the process is completed. This processing is typically carried 
out with a specific algorithm for sequence estimation, for example, the well-known Viterbi 
algorithm. 

[0026] As mentioned above, and as should be apparent from Figs. 3, the performance of 

the equalizer 46 is dependent on the quality of the transmission channel estimation made at 47. 

10 This estimation, in turn, is highly dependent on the accuracy of the burst synchronization and 
equalizer window sizing at 48. The burst synchronization and equalizer window size are 
interdependent on one another, and this interdependency makes efficient optimization not trivial 
in different channel conditions. In contrast to the traditional equaUzer with a fixed span, a 
varying window sizing adaptive to the actual delay spread of the channel is made possible by the 

15 efficient channel estimation proposed here so that a much improved performance can be reached 
in both long delay spread channels such as Hilly Terrain, and short delay spread charmels such as 
Typical Urban. 

[0027] Assuming an equalizer window size m, there are two approaches for burst 

synchronization; a correlation-based approach and a least square error (LSE) approach. 
20 [0028] Synchronization in current GSM receivers is done using a correlation-based 

approach by correlating the central 16 bits of the known 26-bit training sequence with the 
central 16 bits of the training sequence of the received bxirst. 

o 
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20 

c(k)- Sr(i+k)t(iX k=a...,w-l 
i=5 



(1) 



where w is the size of the synchronization window, r(i) is the received samples and t(i) is the 
known training sequence. Synchronization can be determined by either center of gravities of or 
maximimi energy in the equahzer span of c(k). 
5 [0029] hi the LSE approach, for each possible synchronization point, a least square 

estimation of the channel is made, 

H = ^-'T^R = {T^T)-' T^R (2) 

where R is the vector of received symbols, and O is an mxm auto-correlation matrix of the 
training sequence symbols, 
10 <D = r^r (3) 

T is a Toeplitz matrix of the training symbols 



t t , 

m m-\ 



(4) 



1 5 Synchronization is determined by least square error criteria 



p = argmm 



V 



(5) 



where a> \ is a penalty factor, {k) is an estimate of the received symbols 

m-I 



/=0 



(6) 



[0030] Neither of the above approaches, however, is fully satisfactory. The correlation- 

20 based algorithm suffers from performance degradation due to inaccurate synchronization. 
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Performance losses are especially significant in long dispersive channels. The LSE-based 
algorithm provides higher quality channel estimation when accurate synchronization is obtained 
with equation (5); however, it has not been widely implemented due to high computational 
complexity, since, in general, equation (6) requires knowledge of all the channel taps at every 
possible synchronization position with sufficient accuracy. The complexity is even higher if 
channel estimation has to be adapted for different equalizer window size m using equation (2). 
[0031] In a traditional GSM receiver, only 5-tap channel estimation is usually made, 

regardless of the actual channel conditions. It is known that if the 16-symbol mid-amble of the 
training sequence is used for a 5-tap channel estimation, 



10 



r(5,i) = 



^20 ^19 



(7) 



the auto-correlation matrix of the training sequences have a special form 



0(5,l) = r^ (5,1)7(5,1) =16/5 (8) 

15 

where h is a (5 x 5) identity matrix. Here we have made the channel-taps (the order of channel 
vector polynomial), which determine the row dimension of and the index of the first training 
symbol in T explicit parameters of the function. Equation (8) has been used in 5-tap channel 
estimation, 

20 H = 0-'T^R--—T^R (9) 

16 
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In general, however, when Equation (2) is used for initial channel estimation for an equaUzer 
window other than 5, matrix inversions is required. 

[0032] We have now discovered that the GSM/EDGE training sequences have two 

additional properties, not previously exploited, which permit channel estimation to be carried out 
5 more efficiently. More particularly, rewriting equation (8) in the following more general form 
reveals that the autocorrelation matrix of any consecutive 16-symbol segment of the GSM/EDGE 
training sequence is, within certain range, both shift invariant and order invariant. 
Property I. Shift Invariant: 

0(m,d)^T^(m,d)Tim,d) = l6[^ 0<d<{26-i\6+m)+\) (lO) 

10 Property II, Order Invariant: 

im,d)=T^{m,d)T{m,d)=l6l^ l<m<7 (11) 

Property I is due to the cycUc prefix nature of GSM/EDGE training sequences, in which the 26- 
15 bit training sequence is constructed by prepending the last 10 bits to its front as shown in Fig. 2, 
This property enables a channel estimation of the form of equation (9) be carried out on delayed 
(shifted) training sequence segments. This flexible delay brings two important benefits: 

(1) ISI corrupted leading symbols can be avoided in the computation; and 

(2) the leading taps of a channel can be estimated using the same training sequence 
20 segments, regardless of the size of the equalizer window. 

[0033] Property II enables channel estimation to be done for different equalizer windows 
ranging from 1 to 7 taps. It can be stretched fiirther. It has been found that when /w = 8, four of 
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the eight specified training sequence code (TSC) satisfy equation (11), and the auto-correlation 
matrix for the other four training sequences is in a specific sparse form: 

(12) 



m,d)= 



16 0 
0 16 



0 ■ ■ 16 ■ 
±4 0- -16 



The inverse of equation (12) has the same sparse, symmetric structure: 



16 



a 0 
0 1 



a 



(13) 



where a = 1 .0667, p = ± 0.2667. This means that if equation (9) is used for an estimation of 8- 
tap channels, only minor justification is necessary for certain TSCs: 



'K' 




a 




'K 


A. 






a 





(14) 



[0034] Property II indicates that not only a 5 -tap channel, but all the 1-8 tap channels can be 
estimated efficiently using equation (9) without matrix inversion if any consecutive 16-symbol 
segment of the training sequence is used. Utilization of these properties has a fundamental 
impact on the computational complexity of the channel estimation for different equalizer window 
size. 

[0035] When an LSE approach is used for joint synchronization and equalizer window 
sizing, for a hypothetical sync position p and an equalizer window size m, by utilizing the 
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properties discussed above, channel estimation can be computed as in equation (9). With the 

exception of the constant factor, it can be further expressed in an explicit convolution form, 

4^ 0<p<w 
hip,k) = yt(m-\-i-k)rim + i-^p), (15) 

0<k <m 

This can be considered as a FIR-style computation. For the next tap with A:+l, a recurrent 
5 relation can be easily established, 

15 

h(p,k + 1) - Z t(m+i-(k+l))r(p+m4-i) 

(16) 

_ ^t(m+i-k)r(p4-l+m+i) 
i=-l 

= li(p+ l,k)+ t(m- 1- k)r(pH- m) - t(m+ 15- k)r(p + m+ 16) 

Using this recursion in an IIR-style computation, a new tap can be calculated using only 4 real 
MAC operations (since the training sequences are real). This two-dimensional (sync-point tap- 
10 position) recursive relation can be depicted as shown in Fig. 5 for m = 5, w = 8. Computation of 
the last row and the last column is done by using equation (15), the rest can be recursively 
calculated using equation (16). 

[0036] To avoid channel estimation error caused by ISI corrupted receiving training 
symbols, by property I, we can always choose a training sequence segment that is close to the 
15 end. In practice, we often choose the last but one m + 15 symbols of the training sequence in the 
estimation (the final symbol in the training sequence is not used, to avoid possible precursor ISI 
corruption). When the next last m+l5 symbols are chosen, equation (15) can be rewritten in a 
form independent of the equalizer window size m 

h(p,k)^ I t(i-k^9)r(pU+9), oJf<M on ^^^"^ 

i - 0 rnax 

Equation (16) can be adjusted accordingly by simply replacing m with 9. 
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[0037] With the present invention, an efficient initial channel estimation method has been 
obtained for different equalizer window size spanning from 1-8 taps by utilizing the unique 
properties of GSM/EDGE training sequences described above. An efficient channel estimation 
at different sync points for different window size is essential to a practical implementation of 
5 joint optimization of synchronization and equalizer wmdow sizing. The computational 
complexity of initial channel estimation by the approach according to the present invention is 
less than 10% of that of a brute force LSE approach. 

[0038] The presented method has been put into simulation and verified. Simulation shows 
that a significant performance gain is observed in long time dispersive channels with joint 

1 0 optimization compared to the correlation based synchronization. 

[0039] It should be emphasized that the term comprises/comprising when used in this 
specification is taken to specify the presence of stated features, integers, steps or components but 
does not preclude the presence or addition of one or more other features, integers, steps, 
components or groups thereof. 

15 [0040] It should also be emphasized that while what has been described herein constitutes 
presently preferred embodiments of the invention, it should be recognized that the invention 
could take numerous other forms. Accordingly, it should be understood that the invention should 
be limited only insofar as is required by the scope of the following claims. 
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